Role of the conduction electrons in mediating exchange interactions in Heusler alloys 
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Because of large spatial separation of the Mn atoms in Heusler alloys {dMn-Mn > 4 A) the Mn 
3d states belonging to different atoms do not overlap considerably. Therefore an indirect exchange 
interaction between Mn atoms should play a crucial role in the ferromagnetism of the systems. To 
study the nature of the ferromagnetism of various Mn-based semi- and fuU-Heusler alloys we per- 
form a systematic first-principles calculation of the exchange interactions in these materials. The 
calculation of the exchange parameters is based on the frozen-magnon approach. The Curie temper- 
ature is estimated within the mean-field approximation. The calculations show that the magnetism 
of the Mn-based Heusler alloys depends strongly on the number of conduction sp electrons, their 
spin polarization and the position of the unoccupied Mn 3d states with respect to the Fermi level. 
Various magnetic phases are obtained depending on the combination of these characteristics. The 
magnetic phase diagram is determined at zero temperature. The results of the calculations are in 
good agreement with available experimental data. The Anderson's s-d model is used to perform a 
qualitative analysis of the obtained results. The conditions leading to diverse magnetic behavior are 
identified. If the spin polarization of the conduction electrons at the Fermi energy is large and the 
unoccupied Mn 3d states lie well above the Fermi level, an RKKY-type ferromagnetic interaction is 
dominating. On the other hand, the contribution of the antiferromagnetic superexchange becomes 
important if unoccupied Mn 3d states lie close to the Fermi energy. The resulting magnetic behav- 
ior depends on the competition of these two exchange mechanisms. The calculational results are in 
good correlation with the conclusions made on the basis of the Anderson s-d model which provides 
useful framework for the analysis of the results of first-principles calculations and helps to formulate 
the conditions for high Curie temperature. 

PACS numbers: 75.50.Cc, 75.30.Et, 71. 15. Mb 



I. INTRODUCTION 

In recent years, Heusler alloys have become the sub- 
ject of intensive experimental and theoretical investiga- 
tions. The strong interest to these systems is mainly 
due to two unique properties: half-metallic behav- 
ior and martensitic phase transformationSfi"- The half- 
metallicity was first predicted by de Groot and collab- 
orators in 1983 when studying the band structure of a 
half- Heusler alloy NiMnSb."^ Then the half-metallic ferro- 
magnets have become one of the most studied classes of 
materials^i^i^ii2i9iioaiii2ii3^i4j^ rpj^^ existence of a gap 

in the minority-spin band structure leads to the 100% 
spin polarization of the electron states at the Fermi level 
and makes these systems attractive for applications in the 
emerging field of spintronics. Besides strong spin polar- 
ization of the charge carriers, the half-metallic materials 
should have a crystal structure compatible with the in- 
dustrially used zincblende semiconductors and possess a 
high Curie temperature to allow the applications in the 
devices operating at room temperature. The available 
experimental information shows that Heusler alloys are 
promising materials also in this respecti ^^i^^i^^i^^ 

At low temperatures several Heusler compounds (i.e., 
Ni2MnGa, Co2NbSn) undergo a structural transforma- 
tion from a highly symmetric cubic austenitic phase to a 
low symmetry martensitic phase. The compounds that 
are magnetic in the martensitic phase can exhibit two 
unusual effects: magnetic shape memory (MSM) effect 



and inverse magneto-caloric effect i^i^i In MSM alloys 
an external magnetic field can induce large strains when 
applied in the martensitic state. The MSM alloys are 
of great interest as promising smart materials for future 
technological applications i^°i^'^i^'^'^^ They can be used as 
sensors and actuators in different fields of applications. 
The inverse magnetocaloric effect (IMCE) has its origin 
in martensitic phase transformation that modifies the 
exchange interactions due to the change in the lattice 
characteristics. The inverse MCE effect was reported for 
samples with compositions close to Ni2MnZ (Z = Ga, 
Sn) i^iS£i^i2^i22, In the martensitic phase, an adiabatic 
application of a magnetic field, rather than removal of 
the field as in ordinary MCE, causes the sample to cool. 
This feature is regarded promising for the development of 
economical and ecological refrigerants working near room 
temperature as an alternative to conventional vapor-cycle 
refrigeration. 

Besides being promising materials for various appli- 
cations the Heusler alloys constitute a class of systems 
that is important also for fundamental researches. A 
wide diversity of magnetic properties makes Heusler al- 
loys critical test systems for the theoretical models of ex- 
change interactions. Indeed, within the same family of al- 
loys one finds very different magnetic behavior: itinerant 
and localized magnetism, ferrimagnetism, antiferromag- 
netism, helimagnetism and other types of non-coUinear 

ordering 30-31-32.33.34.35.36 

Despite a key role of the magnetism in the properties 
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of Heusler alloys experimental and theoretical studies of 
the exchange interactions in Heusler alloys are still rare. 
The first important information on the exchange cou- 
pling in these systems was provided by the inelastic neu- 
tron scattering experiments of Noda and Ishikawa, and 
Tajima et. ai, in the late seventies^i^. The authors 
measured the spin wave spectra of Ni2MnSn, Pd2MnSn 
and Cu2MnAl for various directions in the Brillouin zone 
and analyzed the results of the measurements within the 
Heisenberg model. They obtained a long-ranged and os- 
cillatory behavior of the exchange interactions. The os- 
cillations were reaching beyond the eight nearest neigh- 
bors in all three compounds. This behavior of the ex- 
change interactions was considered as an evidence for 
an indirect exchange coupling, mediated via conduction 
electrons. The results were interpreted using either an 
RKKY model or s-d mixing model of Anderson (double 
resonance model) j^i^SdldS The double resonance model 
was found to be more relevant for Heusler alloys due to 
strong mixing of the Mn 3d states with the conduction 
electron states of non magnetic 3d and sp atoms. 

This initial attempt of the theoretical analysis did not 
succeed in the description of the sign and the magnitude 
of exchange interactions between nearest and next near- 
est neighbors in contrast to good agreement for larger 
distances4^ The failure of the theories to describe near- 
neighbor exchange interactions was attributed to the 
asymptotic approximations. Price has shown that An- 
derson s-d mixing model free of asymptotic approxima- 
tions was able to capture qualitative features of the ob- 
served spin wave spectra of Pd2MnSn4i Malmstrom and 
Geldart discussed the effect of the finite spin distribu- 
tion around Mn atoms on the RKKY interactions."*^ The 
authors showed that, in spite of simplified treatment of 
electron band structure, finite spin distribution could 
bring the calculated values of the exchange interactions 
in agreement with the exchange interactions determined 
for the Ni2MnSn and Pd2MnSn compounds from the ex- 
periments. 

An important feature of the model-Hamiltonian ap- 
proaches is the possibility of separate study of differ- 
ent exchange mechanisms. However, the use of ad- 
justable parameters and strong simplification of the 
band-structure strongly restrict the possibility of reliable 
predictions for concrete materials. The development of 
the parameter-free density functional theory (DPT) has 
played a crucial role in the understanding of the physical 
properties of itinerant-electron ferromagnets. 

The first contribution to the density functional theory 
of the exchange interactions and Curie temperature in 
Heusler alloys was made in an early paper by Kiibler and 
collaborators where the microscopic mechanisms of the 
magnetism of these systems were discussed on the basis 
of the comparison of the ferromagnetic and antiferromag- 
netic configurations of the Mn moments.— By analyzing 
the Mn 3 c! density of states for different magnetic config- 
urations the authors proposed the mechanism of an indi- 
rect exchange coupling between Mn moments. Recently, 



the studies of the interatomic exchange interactions and 
Curie temperatures in Heusler compounds were reported 
by the present authors and Kurtulus et a;, 47^48,49^50,5i 



II. AIMS AND STRUCTURE OF THE PAPER 

Our previous studies on experimentally well estab- 
lished Ni-based compounds Ni2MnZ (Z = Ga, In, Sn, Sb) 
revealed a complex character of the magnetism in these 
systems4i In particular, the obtained long range and os- 
cillatory behavior of the exchange interactions as well 
as their strong dependence on the sp atom (Z) gave an 
evidence for the conduction-electron mediated exchange 
interactions in Heusler alloys. The importance of the sp- 
electrons in the formation of the magnetic properties has 
been demonstrated by first-principles calculations also 
for a number of other systems.— i^*^ The experimen- 
tal studies as well pointed out to the important role 
of the sp-electrons. The early measurements by Web- 
ster et al. have shown this for quaternary Heusler alloys 
Pd2MnIni_a;Sna; and Pd2MnSni_2;Sba;.— Recent experi- 
mental studies on the Mn-based semi Heusler compounds 
Nii_a;Cua;MnSb and AuMnSni-jjSb^: revealed similari- 
ties in the behaviour of the two systems with the variation 
of the number of valence electronsi^i^ In particular, the 
Curie temperatures of both systems decrease by the same 
amount with increasing concentration x and both sys- 
tems have similar Tc values for compositions with equal 
numbers of valance electrons. These studies have shown 
that both the non-magnetic 3(i-atoms (X) and the sp- 
atoms (Z) play important role in establishing magnetic 
properties. The explanation of the observed trends in 
magnetic characteristics of semi-Heusler alloys requires 
consideration of competing exchange mechanisms con- 
tributing to the formation of the magnetic state. One 
of the important aims of this work is to show that the 
results of the parameter-free DPT calculations for many 
systems can be qualitatively interpreted in terms of the 
competition of two exchange mechanisms. 

In contrast to our previous studies, here we go beyond 
the stoichiometric compositions and treat several Mn- 
based semi and full Heusler alloys within virtual crystal 
approximation in order to understand the dependence of 
physical characteristics on the valance electron number. 

The qualitative interpretation of the calculational re- 
sults is based on the analysis of the Anderson's s-d mix- 
ing model42. Since the atomic Mn moments are well de- 
fined the application of the model is well founded. Por 
the last four decades the Anderson model has been suc- 
cessfully applied to variety of problems of condensed 
matter physics and contributed to better understand- 
ing of numerous experimental properties. The systems 
studied within the Anderson model vary from dilute 
3d magnetic impurities in non-magnetic metals^'^'^'^ 
to complex systems like Heusler alloys^^iiiii^, rare 
earths^- magnetic multilayer a^^^i^^i^^ and diluted mag- 
netic semiconductors i^i^iSSiilS 
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FIG. 1: Clb and L2i structures adapted by the half- and fuU- 
Heusler alloys. The lattice consists from 4 interpenetrating fee 
lattices. In the case of the half-Heusler alloys (XYZ) one of the 
four sublattices is vacant. In VGA, the Z site is occupied by 
a pseudoatom with a fractional number of valence electrons. 



TABLE I: Experimental and theoretical lattice parameters in 
XkMnZ ( X = Pd, Gu; k = 1, 2; Z = In, Sn, Sb, Te). For 
explanations see text. 



In the present paper we show that the diversity of 
magnetic behavior in Mn-based Heusler alloys can be 
interpreted in terms of the competition between ferro- 
magnetic RKKY-type exchange and antiferromagnetic 
superexchange. A special attention is devoted to the role 
of the sp and non-magnetic 3d atoms in mediating ex- 
change interactions between Mn atoms. The spin polar- 
ization of the conduction electrons appears to be one of 
the key parameters in the formation of magnetic char- 
acteristics. We obtain strong correlation between spin 
polarization of the sp electrons and the strength of the 
exchange interactions and, as a result, the value of the 
Curie temperature. It is shown that the position of un- 
occupied peaks of the Mn 3d density of states plays an 
important role in the determining the value of the anti- 
ferromagnetic superexchange coupling while the proper- 
ties of the conduction-electron states at the Fermi level 
are mainly responsible for the ferromagnetic RKKY-type 
exchange interaction. We compare the influence of non- 
magnetic 3d versus Ad atoms on the exchange coupling 
by considering the stoichiometric Ni2MnZ and Pd2MnZ 
(Z=In, Sn, Sb, Te). 

One of the main conclusions of the present treatment 
is that numerous features of the magnetism of the Mn- 
based semi and full Heusler alloys with different chem- 
ical composition can be described with few parameters. 
Therefore, the tuning of these parameters can be con- 
sidered as a tool in the fabrication of the materials with 
desired magnetic properties. 

A part of the present work has been published 
elsewherci^ The remaining of the paper is organized as 
follows. In Section III we present the details of the calcu- 
lational approach. Section IV contains the discussion of 
the density of states and magnetic moments. In section 
V we dwell on the exchange coupling mechanisms and 
Section VI gives the conclusions. 
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III. TECHNICAL DETAILS 

Semi- and fuU-Heuslcr alloys crystallize in the Clb and 
L2i structures respectively (Fig. [T]). The lattice consists 
from 4 interpenetrating fee lattices. In the case of the 
semi-Heusler alloys (XYZ) one of the four sublattices is 
vacant. The Bravais lattice is in both cases fee. 

The calculations are carried out with the aug- 
mented spherical wave method (ASW)ii within the 
atomic-sphere approximation (ASA).''^ The exchange- 
correlation potential is chosen in the generalized gradient 
approximationjS A dense Brillouin zone (BZ) sampling 
30 X 30 X 30 is used. The radii of all atomic spheres are 
chosen equal. In the case of semi-Heusler alloys we intro- 
duce an empty sphere located at the unoccupied site. 

We focuss on the systems where the total magnetic 
moment is confined to Mn sublattice that simplifies the 
interpretation of the obtained results. To this end we 
consider Pd and Cu containing Mn-based semi and full 
Heusler alloys that can be written in a compact form as 
follows 



XkMnZi_mZ' 



< m < 1, 
< m < 1, 



k = 1 for s-H 
k = 2 for f-H 



where X = (Pd,Cu), (Z, Z') = (In,Sn), (Sn, Sb), (Sb, 
Te) and s-H (f-H) stands for semi-Heusler (fuU-Heusler) 
alloys. To account for non-integer electron numbers we 
use the virtual crystal approximation (VCA).^"' In VGA, 
the Z site occupied by In and Sn atoms according to their 
concentration which is described by an atom with frac- 
tional number of electrons (1 — m) z^" + mz^" where 
z"^" is the number of electrons of Sn and similarly for 
other chemical elements. In this approximation only the 
spin magnetic moment and the density of states (DOS) 
of the pseudoatom can be calculated. The properties of 
the pseudoatom cannot be projected on the constitut- 
ing atoms. An advantage of VGA is the possibility of 
a continuous variation of the electron number without 
resorting to large super cells. 

In the last column of Table U we present the lat- 
tice parameters used in the calculations. The remaining 
columns give the lattice constants of the experimentally 
existing systems. Some of the compounds are not yet 
synthesized. For the synthesized systems, the alloys con- 
taining the sp-atom (Z) from the same row of the Periodic 
Table have similar values of lattice constants. Therefore, 
we use the average lattice parameters of the synthesized 
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systems in the calculations for all Z constituents. As seen 
from Table [I] the difference between experimental lattice 
parameters and the parameters used in the calculations 
is less than 0.5%. 

To determine the interatomic exchange interactions we 
use the frozen-magnon technique^ and map the results of 
the calculation of the total energy of the helical magnetic 
configurations ^^'^^ 

s„ — [cos(qR„) sm9, sin(qR„) sinO, cos9] (1) 

onto a classical Heisenberg Hamiltonian 

Heff ^ -^JijSiSj (2) 

where Jy is an exchange interaction between two Mn 
sites and is the unit vector pointing in the direction of 
the magnetic moment at site i , R„ are the lattice vectors, 
q is the wave vector of the helix, 9 polar angle giving 
the deviation of the moments from the z axis. Within 
the Heisenberg model the energy of frozen-magnon 
configurations can be represented in the form 

E{9,q)=Eo{9)-sinHj{q) (3) 

where Eq does not depend on q and J(q) is the the 
Fourier transform of the parameters of interatomic ex- 
change interactions: 

•^(q) = "^OR- exp(iq • R). (4) 

R 

Calculating E{9, q) for a regular q-mesh in the Bril- 
louin zone of the crystal and performing back Fourier 
transformation one gets exchange parameters Jqr be- 
tween pairs of Mn atoms. 

The Curie temperature is estimated within the mean- 
field approximation (MFA) 

kBTS'^^ = lY.Jo, (5) 



IV. DENSITY OF STATES AND MAGNETIC 
MOMENTS 

The electronic and magnetic structures of Hcusler al- 
loys have been extensively studied earlier and the reader 
is referred to Ref. and the references therein for a de- 
tailed overview. Here we present a brief description of 
the calculational results aiming to provide the basis for 
the discussion of exchange mechanisms in subsequent sec- 
tions and to allow the comparison with previous work. 

A. Density of states 

In this section, we discuss the density of states for 
the stoichiometric compositions of both families of com- 
pounds [XkMnZ ( X = Pd, Cu; k = 1,2; Z = In, Sn, 



Sb, Te)]. The results are presented in Fig. H The DOS 
for non-stoichiometric compositions assumes intermedi- 
ate values. In agreement with the commonly accepted 
picture of the magnetism of the Mn-based Heusler alloys 
we obtain a strong localization of the magnetization on 
the Mn sublattice with a value of the Mn moment close 
to 4/iB. In the following we will show that the value 
of the Mn moment is very robust with respect to the 
variation of the magnetic structure. Therefore we can 
conclude that Mn-based Heusler alloys possess a well- 
defined atomic Mn moment. The robust character of the 
Mn moment results from the large exchange splitting of 
the Mn 3d states. Important that the Mn 3d states of 
only one spin projection (spin-up) are strongly occupied. 
The main part of the spin-down Mn 3d states lies above 
the Fermi level. 

The analysis of the DOS of different compounds reveals 
a relative shift of the Fermi level to a higher energy po- 
sition in the sequence In-Sn-Sb-Te (Fig. [J) that results 
from the increasing number of valence electrons within 
this series. Indeed, a Te atom has three more valence 
electrons than In, two more electrons than Sn and one 
more electron than Sb. The Mn DOS for semi-Heusler 
compounds shows an interesting property that different 
compositions give similar features if the total number of 
the sp-electrons coming from different atoms is the same. 
This similarity is more pronounced near the Fermi level. 
For example, the peaks in the Mn spin down states of 
PdMnSn and CuMnIn are very similar to each other. As 
a consequence, the Mn atomic moments are also simi- 
lar. This property can be related to the symmetry of the 
wave functions in Clf,-type crystal structure of the semi- 
Heusler alloys.'* In the following sections we show that the 
position of the occupied Mn 3 c! states with respect to the 
Fermi level, the sp-electron DOS at the Fermi energy, and 
the structure of the unoccupied part of the DOS are im- 
portant for the understanding of the exchange coupling 
mechanisms in Heusler alloys. 



B. Magnetic moments and local moment behavior 

Most of the Mn-based systems are believed to possess 
well-defined local atomic Mn moments which are usually 
close to AfiB and do not change substantially when go- 
ing from the ordered phase to the paramagnetic state. 
The available inelastic neutron scattering experiments 
and magnetization measurements support this point of 
view. Early studies of the paramagnetic phase of sev- 
eral Mn-based compounds have established the absence 
of spatial magnetic correlations (spin waves), and have 
shown the value of the atomic moment to be in agree- 
ment with the moment obtained from the static suscep- 
tibility measurements Recently, Plogmannet. al., gave 
a detailed study on the degree of magnetic moment lo- 
calization in various Mn-based fuU-Heusler alloys using 
x-ray photoelectron spectroscopy (XPS) and x-ray emis- 
sion spectroscopy (XES) techniques. '^^ Depending on the 
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PdMnZ CuMnZ Pd^MnZ Cu^MnZ 




E(Ry) E(Ry) E(Ry) E(Ry) 

FIG. 2: (Color online) Left panel: Spin-projected atom-resolved density of states of PdMnZ and CuMnZ (Z=In, Sn, Sb, Te) 
for stoichiometric compositions. The shadow areas show the DOS of the Z constituent. The broken vertical lines denote the 
Fermi level. Right panel: The same for full Heusler compounds Pd2MnZ and Cu2MnZ. 



atomic number of the Z element, the authors obtained 
from XES experiments an increasing localization of the 
Mn-3rf states that is related to the larger interatomic 
distances for heavier sp (Z) elements. 

On the theoretical side, in 1984 Kiibler et ai, gave 
a detailed account of the formation of local moments in 
various Mn-based full Heusler alloys using first-principles 
calculations!^ Since then many authors have studied var- 
ious semi- and fuU-Heusler compounds and came to sim- 
ilar conclusion. An interesting features arising from the 
calculations of Heusler alloys is the resemblance of a num- 
ber of physical characteristics such as the position and 
the width of the Mn 3d peaks, exchange splitting and 
the value of the Mn magnetic moment to the correspond- 
ing quantities for diluted Mn impurities in non-magnetic 
metals.~ This observation gives further evidence for the 
localized nature of the magnetism moments in these sys- 
tems. 

In Fig. [3] we present calculated atom-resolved and to- 
tal magnetic moments for both families of Heusler alloys 
PdMnZ, CuMnZ, Pd2MnZ and CusMnZ as a function of 
the sp-electron number of the Z constituent. As men- 



tioned above, the magnetic moment is mostly confined 
to the Mn sublattice. Small moments are induced on 
the Cu and Pd atoms. These induced moments are pos- 
itive (that is parallel to the Mn moments) in a broad 
composition interval. The induced moment of the Z ele- 
ment is negative. The induced moments as a function of 
the Z constituent follow closely the behavior of the Mn 
moment. This correlation is especially well pronounced 
for semi-Heusler alloys (Fig. [3]). The comparison of the 
data for PdMnZ and CuMnZ shows that for equal total 
number of valence electrons the induced moments have 
similar values. In both semi-Heusler systems, the vari- 
ation of the total magnetic moment with Z constituent 
is large. In PdMnZ, it is about 1 fiB- The half of the 
variation comes from the change in the Mn moment. In 
full Heusler alloys the situation is different. Neither a 
correlation of the characteristics of the compounds with 
the same number of the valance electrons nor a substan- 
tial change in magnetic moments with variation of Z is 
obtained. 

Since an increase of the temperature leads to increasing 
deviation of the atomic moments from the direction of the 



6 



PdMnZ 



CuMnZ 




Pd^MnZ 



Cu^MnZ 



0.16 

0.08 

2" 

I -0.08 
o 

s 

u 

•■S 4.4 
e 

4 
3.8 



□ □ Pd 
■ ■ Z 



O-O Mn 
• • Cell 



□ ■ □ Cu 

■ ■ z 



O-O Mn 
• • Cell 



(b) 



In 



Sn 



Sb Te In Sn 



Sb Te 



FIG. 3: (Color online) (a) Calculated atom-resolved and total spin moments (in /xb) in PdMnZ and CuMnZ as a function of 
the sp-electron number of the Z constituent, (b) The same for full Heusler alloys Pd2MnZ and Cu2MnZ. 
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FIG. 4: (Color online) (a) Upper panel: Calculated Mn spin moment in PdMnSn and CuMnSn as a function of 9 for spin 
spiral with q = (0 |), (0 1) in units of 2-k ja. Lower panel: The corresponding total energies AE{0, q) = E{6, q) — E{0, 0). 
For comparison, the results of the force theorem calculations (broken lines) are presented. Jq(S) stand for Jq x M|, M is the 
magnetic moment, (b) The same for Pd2MnSn and Cu2MnSn full Heusler compounds. 



net magnetization it is important to study the properties 
of the noncollinear magnetic configurations. In Fig.lH we 
present the results of such calculations for four stoichio- 
metric compounds: PdMnSn, CuMnSn, Pd2MnSn and 
Cu2MnSn. The results are shown for the spiral struc- 
tures with two different wave vectors q = (0 i) and 
q = (0 1) in units of 27r/a and the polar vector 9 vary- 



ing in the interval from to 90° (sec Eq. [T]). 

The upper panels of Fig. 2] show the 9 dependence of 
the Mn moment whereas the bottom panels present the 
9 dependence of the total energy. Note that by the vari- 
ation of the wave vector of the spiral and of angle 9 the 
magnetic structure can be continuously transformed from 
ferromagnetic to antiferromagnetic. The spirals with 
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9 = and arbitrary q correspond to the ferromagnetic 
state whereas the structure with 6 = 90° and q = (0 01) 
is antiferromagnetic. 

The analysis of Fig. 2] shows that for all four com- 
pounds the magnetic moment of the Mn atom is prac- 
tically insensitive to the magnetic configuration. The 
relative change of the Mn moment in transition from fer- 
romagnetic to antiferromagnetic state is less than 3%. 
Therefore the treatment of the Mn moment in these sys- 
tems as a local property of the Mn atom is well founded. 

The energy of the spiral structures increases 
monotonously with increasing 9 and reaches the maxi- 
mal value for the antiferromagnetic state. The minimum 
at 6* = shows that the ground state is ferromagnetic. 
In Fig. m we compare the total energy calculated self- 
consistently with an approximation to the total energy 
obtained with the application of so-called magnetic force 
theorem^® The agreement between two types of calcula- 
tions is good in the whole 9 interval that is characteristic 
for systems with well defined atomic moments and allows 
the use of the force theorem for magnetic configurations 
deviating strongly from the ground state. 



V. EXCHANGE COUPLING MECHANISM 

In this section we suggest an interpretation of the re- 
sults of the first-principles calculations on the basis of 
the Anderson s-d mixing model. The discussion is di- 
vided into four parts. In the first part, the calculated 
Heisenberg exchange parameters are presented. In the 
second part we briefly describe the Anderson s-d mixing 
model and the exchange mechanisms resulting from the 
perturbative treatment of this model. The role of the 
conduction-electron spin polarization in exchange cou- 
pling as well as the contribution of the sp and non- 
magnetic 3d atoms to the formation of the magnetic state 
is discussed in the third part. The last part is devoted to 
the study of the influence of non-magnetic 3d versus Ad 
atoms on the exchange coupling. 

A. Heisenberg exchange parameters 

In Fig. O we present calculated Mn-Mn exchange pa- 
rameters for six nearest neighbors as a function of sp- 
electron concentration. In the upper-right corner of each 
panel the number of atoms in the corresponding coor- 
dination sphere is given. The exchange parameters for 
larger interatomic distances 10 a) are shown in Figs.lH] 
for selected compounds PdMnIn, CuMnIn, Pd2MnIn and 
Cu2MnIn. The absolute value of the exchange parame- 
ters decays quickly with increasing interatomic distance 
and the main contribution to Tc comes from the interac- 
tion between atoms lying closer than 3a. No sizable con- 
tribution is detected after 5a. However, the RKKY-type 
oscillations become visible up to very large interatomic 
distances when the exchange parameters are multiplied 



by (d/a)^ (Fig- E]) where d is the distance between the 
interacting Mn atoms and a is the lattice constant). 

In agreement with the results of our previous calcula- 
tions on Heusler alloys a strong dependence of the ex- 
change parameters on the Z constituent for both families 
of systems is obtained. As seen from Fig. [S]all exchange 
parameters oscillate between ferromagnetic and antifer- 
romagnetic values with increasing sp-electron concentra- 
tion. As discussed below these oscillations are related to 
the properties of the electron structure of the systems. 

Considering two nearest neighbor exchange parame- 
ters we notice that they have ferromagnetic character 
for a broad composition range and dominate over the 
rest of parameters. The remaining parameters are much 
weaker. The first two nearest neighbor exchange interac- 
tions are responsible for very high Curie temperatures in 
Cu-based full Heusler alloys as well as in both classes of 
semi Heusler alloys. 

A distinct feature of the exchange interactions in semi- 
Heusler alloys is that the maximum of the exchange in- 
teractions for both PdMnZ and CuMnZ corresponds to 
the similar numbers of the sp-electrons. This correlates 
with the conclusions of the preceding sections where a 
similar behavior was obtained for the density of states 
and the magnetic moments. The shift of the maxima for 
two systems is explained by the fact that Pd has one sp- 
electron less than Cu. The properties of the exchange 
interactions are reflected in the properties of the Curie 
temperature (Fig. [7]) where we also obtained a relative 
shift of the maxima of the two curves corresponding to 
one sp-electron. However, no such correlation is obtained 
for the full Heusler compounds. 



B. Indirect exchange coupling: RKKY-type 
exchange and Superexchange 

In previous calculations on Co2MnSi and Mn2VAl we 
have shown that in the systems with several magnetic 
sublattices the direct exchange coupling Mn-Co or Mn-V 
between neighboring 3d atoms can dominate over the in- 
direct Mn-Mn interactions."'^ For the Mn-based Heusler 
compounds considered in this paper the direct coupling 
does not play substantial role and can be ignored4i 
Therefore, in the following we will consider only the in- 
direct exchange coupling. 

The DFT is not based on a model Hamiltonian ap- 
proach and does not use a perturbative treatment. 
Therefore various exchange mechanisms appear in the 
results of calculations in a mixed form that does not al- 
low a straightforward separation of the contributions of 
different mechanisms. In this situation, the model Hamil- 
tonian studies relevant to the problem provide useful in- 
formation for a qualitative interpretation of the DFT re- 
sults. Among such approaches the Anderson s-d mixing 
model is an appropriate tool for Mn-based Heusler alloys 
because of the localized nature of Mn moments in these 
systems. 
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FIG. 5: (Color online) Left panel: First six nearest neighbor Mn-Mn exchange parameters in PdMnZ and CuMnZ as a function 
of the sp-electron number of the Z constituent. Also given are the number of atoms within corresponding coordination spheres. 
Right panel: The same for fuU-Heusler alloys Pd2MnZ and Cu2MnZ. 
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FIG. 6: (Color online) (a) Upper panel: Mn-Mn exchange interactions in PdkMnIn (k — 1, 2) as a function of distance up to the 
10a. Lower panel: RKKY-type oscillations in exchange parameters for corresponding compounds, (b) The same for CukMnIn 
(k=l,2). 



Anderson s-d model for a single 3d impurity embedded Hamiltonians^ 
into a nonmagnetic metallic host is given by the following 



kcr 



'^Vkdia'l^ada + auaa^a) + Uud^Udi 
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where aj^^. (flko-) and a^^ (fldo-) create (annihilate) elec- 
trons with spin a in the band states and on the impu- 
rity, respectively, Uda is the number operator for local- 
ized electrons of spin cr, ?7 is the on-site Coulomb re- 
pulsion between two localized electrons which favors the 
single occupation of the impurity level. Vkd represents 
the coupling between the impurity level and the conduc- 
tion electrons in the metal. This interaction causes the 
mixing of the band states and impurity level. Such a mix- 
ing gives rise to broadening of the impurity states. The 
ratio U /V determines whether a local magnetic moment 
in a metallic host can be formed. The first-principles 
calculations for Sd impurities in non-magnetic metals 
gave large values for the magnetic moments reflecting 
the strong Coulomb repulsion U between Sd orbitals.^^ 
Note that in complex systems like rare earths and Heusler 
alloys one deals with a lattice of magnetic atoms. In 
this case, the relevant model is the periodic Anderson 
s-d modelM>2i^ 

The mixing interaction V induces spin polarization 
in the conduction electron sea and spatial propagation 
of this polarization gives rise to an effective indirect 
exchange coupling between distant magnetic moments. 
In general, this interaction contains two distinct pro- 
cesses: electrostatic Coulomb exchange interaction and 
s-d (or sp-d) mixing interaction. The former induces 
a net positive spin polarization while the contribution 
of the latter is negative and disappears in strong mag- 
netic limit. The indirect exchange coupling has been 
described using various theoretical schemes. The Green's 
function method and the perturbation theory are among 
the most often used techniques. In 1973 da Silva and 
Falicov showed that, in fourth-order perturbation theory, 
the indirect interaction between two magnetic moments 
can be separated into two contributions!^ ^/indirect — 
Jkkky + Js- The first term is of RKKY type and stems 
from the intermediate states which correspond to low en- 
ergy spin excitations of the Fermi sea, that is exitations 
corresponding to electron-hole pair formation with a spin 
flip transition. In reciprocal space this contribution takes 
the following form.^^'Si^es^ 



RKKY 



(q) = 



rii ,712 ,k 



ni A; I 



IK 



n^k' I 



9{ef — enik)0{en2k' ~ £/) 

^n2k' ^nik 



+ C.C. 



where k' = k + q + G and G is the reciprocal lattice 
vector. This coupling is largely influenced by the de- 
nominator (e„2fe' — Snik)- Therefore the topology of the 
Fermi surface is an important factor influencing the form 
of JBKKYi^) with most important contributions com- 
ing from the stationary wave vectors spanning the Fermi 
surface. These stationary wave vectors translate into 
smoothly decaying oscillations in the real space with a 
ferromagnetic bias in the "pre- asymptotic" regime. Also 
the DOS at the Fermi level, the number of conduction 



electrons and their spin polarization play important role 
in determining the strength of this interaction. 

The second term arises from high-energy virtual charge 
excitations in which electrons from local 3d states of the 
magnetic atom are promoted above the Fermi sea. In the 
reciprocal space, the term has the following form that is 
similar to the form of the RKKY-type term 



^s(q) = - J2 



Tl-l ,712 - A; 

Oismk 



iKiifePlKiafe' P 

(e„2fc' - 



£nik - £d 



+ C.C. 



In contrast to the RKKY contribution, Js(q) does not, 
however, depend on the DOS at the Fermi level and the 
topology of the Fermi surface. Since the sum is taken 
over unoccupied states an important role is played by 
the energy position of the unoccupied 3d states of the 
magnetic atom. The closer the states to the Fermi level 
the larger Js (q) ■ This interaction is always antiferromag- 
netic and its strength decays exponentially with distance. 
A broadening of the 3d levels induces weak oscillations 
in this coupling. In addition to these parameters, the 
position of the occupied 3d levels with respect to the 
Fermi energy and the strength of the mixing interaction 
V strongly influence coupling mechanisms. 

As shown in Ref. 66 the q limit simplifies the 
above expressions and is useful for the qualitative analy- 
sis. For q = the RKKY-type term becomes. 

Jrkky(O) = V^DieF)/El 

where D{ep) is the density of states at the Fermi level. 
Eh is the energy required to promote an electron from 
occupied 3d level to the Fermi level. Parameter Eh is not 
well defined in Heusler alloys because of the broadening 
of the Mn 3d levels into the energy bands crossing the 
Fermi level. In all alloys studied the occupied Mn 3d 
peaks lie below —0.1 Ry where the energy is counted off 
the Fermi level (Fig. [2]). On the other hand, the q = 
limit of the superexchange term cannot be expressed in 
terms of the density of states and has more complex form 



nk 



ink — Eh] 



-3 



In addition, we would like to comment on the following 
two points. First, although the perturbative derivation of 
above expressions is based on the assumption of two mag- 
netic impurities embedded into metallic host, the gener- 
alization to the periodic lattices is straightforward. As 
shown by G. Malmstrbm et. at, and Price the final re- 
sults differ by a phase factor..d2jii. Second, there are two 
different limits in the description of the exchange mech- 
anisms within Anderson s-d model. In weak magnetic 
limit where the coupling is dominated by the s-d mixing 
of the local and conduction electron states, both mech- 
anisms mentioned above coexist and their relative con- 
tributions are determined by the details of the electronic 
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structure of the system. In the opposite hmit (strong 
magnetic limit) in which the coupling is primarily due 
to the electrostatic Coulomb exchange interaction, the 
second term is not present and the first term is reduced 
to the conventional RKKY interaction.^^ The systems 
considered in this paper are characterized by strong s-d 
hybridization and therefore we expect the presence of the 
contributions of both exchange mechanisms. 



C. Conduction electron spin polarization 

Now we turn to the interpretation of the DFT results in 
terms of the two exchange mechanisms discussed in the 
previous section. To estimate the relative contribution 
of the ferromagnetic RKKY-type exchange we present in 
Fig. [7] the calculated conduction electron spin polariza- 
tion in both families of Heusler alloys as a function of 
the electron number of the Z constituent. The analysis 
of the calculational data allows us to draw a number of 
important conclusions. First, in both families of alloys 
the direction of the induced spin polarization is opposite 
to the direction of the Mn moment in a broad interval 
of compositions. This feature reveals the primary role 
of the sp-d mixing in the exchange coupling and justifies 
the use of the Anderson s-d model for the description of 
the magnetism in these systems. 

Another remarkable feature is a very clear correlation 
between the spin polarization and the mean-field Curie 
temperature (or exchange parameters) in a large part 
of the phase diagram (Fig. [7]). Indeed, the compounds 
with very large spin polarization are characterized by the 
value of the Curie temperature that is also very high. In- 
terestingly, for the zero polarization the Curie tempera- 
ture also vanishes or assumes very small values refiecting 
the dominating character of the ferromagnetic RKKY- 
type exchange mechanism in establishing magnetic order. 
However, at some regions superexchange mechanism be- 
comes important. This can be seen in Fig. [7] where for 
PdMnIni_a;Sn2. (x < 0.8) system we obtain an antifer- 
romagnetic order in spite of very large spin polarization. 
Further insight into the nature of the coupling can be 
gained from Fig. [2] where we present atom and spin re- 
solved density of states of both families of compounds for 
stoichiometric compositions. As pointed out above, the 
superexchange mechanism is sensitive to the DOS above 
the Fermi energy. As seen in Fig. [2] the Mn 3d states 
provide the main contribution to the DOS in this region. 
For PdMnIn the DOS peak above the Fermi level assumes 
the largest value and, as a result the antiferromagnetic 
superexchange dominates over the ferromagnetic RKKY- 
type exchange giving rise to AFM order. In transition 
from In to Sn the peak gradually decreases and there- 
fore the superexchange becomes less important. Around 
Z = Sb it almost disappears leading to FM order. How- 
ever, another large peak approaches the Fermi level for 
Z = Te that turns the systems from a ferromagnet into 
an antiferromagnet with a non-coUinear ordering in the 



intermediate region due to the competition of two mecha- 
nisms. The situation is very similar in the case of CuMnZ 
where we also obtain a rich magnetic behavior. 

On the other hand, for the fuU-Heusler alloys the mag- 
netic phase diagram is rather simple compared to the 
semi-Heusler systems. In Cu2MnZ, the ground state is 
ferromagnetic for all Z. The mean-field Tc qualitatively 
follows the behavior of conduction electron spin polar- 
ization, assuming largest value for Z = In and gradually 
decreasing in the In-Sn-Sb-Te sequence. Around Sb the 
conduction electron spin polarization increases that in 
principle should lead to the increase of the Curie temper- 
ature according to the discussion of the previous section. 
However, in this region the Tc further decreases reveal- 
ing substantial contribution of the antiferromagnetic su- 
perexchange. This is reflected in position of the Mn 3d 
states that is very close to the Fermi level (Fig. [5]). 

In the Pd-based full Heusler alloys Pd2MnZ the spin 
polarization is rather small, and as a result the Curie 
temperature is low. The spin polarization in In-Sn in- 
terval is independent of the Z constituent, while in the 
rest of the phase diagram it gradually decreases and be- 
comes zero around Te. As for the Tc we obtain a similar 
Z-independent behavior in a large interval of composi- 
tions (from In to Sb). However, from Sb to Te the Tc 
sharply decreases due to dominating contribution of su- 
perexchange. The calculated Tc values are in good agree- 
ment with the experiments for the Pd-based full Heusler 
alloys, while they are overestimated in PdMnSb and Cu- 
based full Heusler compounds. The following explana- 
tion of this property can be suggested. In the Pd-based 
full Heusler alloys the more distant exchange parameters 
contribute with a substantial weight to the formation of 
the Curie temperature (Fig. [5]) whereas in PdMnSb and 
Cu-based full Heusler compounds only the first and sec- 
ond nearest neighbor parameters determine the Tc. In 
the latter case the mean-field approximation is less ex- 
act and overestimates the Curie temperature. For such 
systems, the random-phase approximation (RPA) is ex- 
pected to provide a better description of the Tc. 

The results of the first principles calculations and the 
qualitative analysis on the basis of the Anderson s-d 
model allow us to formulate two conditions for high Curie 
temperature, (i) The conduction electron spin polariza- 
tion should be maximal, (ii) The DOS above Fermi level 
should be minimal. Indeed, as seen from Figs. [5] and [3 
the compounds CuMnSn, PdMnSb and Cu2MnIn that 
satisfy both conditions posses very high Curie tempera- 
tures. The knowledge of the conditions for high Curie 
temperature is an important help in the fabrication of 
the materials with desired properties. From the point of 
view of the second condition the half-metallic ferromag- 
nets have an advantage since the gap in the spin-down 
channel decreases the number of the states just above 
the Fermi level. Indeed, a large number of first-principles 
calculations showed that these materials have very high 
Curie temperatures^i^s^'Sl^^iSO 

As mentioned in preceding part the direction of the 
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FIG. 7; (Color online) (a) Upper panel: The sp-electron spin polarization (in /^b) of the X=Pd,Cu and Z constituents. The total 
polarization is given as the sum of X and Z spin polarizations. Lower panel: Mean-field estimation of the Curie temperature 
in PdMnZ and CuMnZ. For comparison available experimental Tc values taken from RefsT, and[3 are presented. FM, NC and 
AFM stand for ferromagnetic, non-coUinear, and antiferromagnetic ordering, respectively . (b) The same for the fuU-Heusler 
alloys Pd2MnZ and Cu2MnZ. 



conduction electron spin polarization with respect to lo- 
cal moment and its amplitude are important in classi- 
fying ferromagnets. From the obtained results the fuU- 
Heusler alloys can be characterized as weak ferromag- 
nets since spin polarization is negative for all Z con- 
stituents, i.e., sp-d mixing is dominating. As for the semi- 
Heusler alloys the situation is similar to the fuU-Heusler 
alloys except some regions in the magnetic phase diagram 
where spin polarization changes sign: Sb-Te interval in 
CuMnSb and the region around Tc in PdMnZ. A positive 
spin polarization does not directly mean that these alloys 
are strongly ferromagnetic. They are, however, close to 
this limit and the Mn-Mn coupling is primarily due to 
the electrostatic Coulomb exchange interaction. Indeed, 
as seen from Figs. [3] and [7] the alloys with large posi- 
tive spin polarization have large Mn magnetic moments. 
In PdMnTe both quantities assume the largest values: 
niMn — 'iAfiB and mgp = 0.18/iB. The dependences of 
the Mn moment and of the spin polarization on the Z 
constituent are correlated (Figs. [3] and [7]). 

Because of the importance of the conduction-electron 
spin polarization a direct comparison with experimental 
measurement of this quantity is desirable. An impor- 
tant information on conduction-electron spin polariza- 
tion is supplied by the measurements of the hyperfine 
fields at non-magnetic sites (X,Z). The strength of the 
transferred hyperfine fields correlates with the amplitude 
of the s conduction electron spin polarization. Indeed, 
the measurements by Campbell and Khoi et. at, showed 
that maximal s-electron spin polarization is found in the 
systems with high Curie temperatures such as Cu2MnAl 



and Cu2MnIni21i2^ However, the hyperfine fields are sen- 
sitive to only the s electron spin polarization and do not 
give information about the polarization of the p electrons 
that are usually dominating in Heusler alloys. For prob- 
ing the total conduction electron spin polarization the 
magnetic Compton scattering profiles are proved to be 
a useful tool. Using this method, Zukowski et. al. ob- 
tained recently a large sp-electron spin polarization in 
Cu2MnAl which is antiferromagnetically coupled to the 
Mn momentsi^ A similar result is obtained by Deb et. 
al., for Ni2MnSnii2i Our calculations are in agreement 
with both experiments. 



D. 3d versus 4d electrons 

In majority of the Heusler alloys with chemical formula 
X2YZ and XYZ the Y site is occupied by the Mn atom, 
while for the X site there is much more freedom: here 
can be any element from 3d, 4d, or bd late transition 
metals (i.e. Fe, Ru, Co, Rh, Ni, Pd, Pt, Cu, Ag and 
Au) . The magnetic moment of the X atom can for many 
alloys be neglected. The exceptions are some of the 3d 
atoms, e.g., Fe or Co. Since the delocalization of the d 
states increases with transition from 3d elements to 4d 
and further to 5d elements one can expect a substantial 
dependence of the magnetic properties of the Heusler al- 
loys for the occupation of the X site with atoms from 
different d series. The purpose of this section is to study 
this dependence by means of comparison of two systems: 
Ni2MnZ and Pd2MnZ (Z = In, Sn, Sb, Te). 
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FIG. 8: (Color online) Left panel: The parameters of the Mn-Mn exchange interactions in Ni2MnZ and Pd2MnZ (Z=In, Sn, 
Sb, Te) for two different lattice constants. Right panel: Spin projected total DOS (left hand side) and sp-electron DOS of X 
and Z atoms (right hand side) for Ni2MnZ and Pd2MnZ (Z=In, Sn, Sb, Te) for the lattice constants of 6.38A. 



TABLE II: Magnetic moments, sp-electron polarization (in fis) and Curie temperatures in Ni2MnZ and Pd2MnZ (Z = In, Sn, 
Sb, Te) for two different lattice parameters. Negative values of Tc mean that ground state is antiferromagnetic. 
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T^"'*(K) 
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Mn 
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Cell 


Tg^*'*(K) 


Ni2MnIn 


0.27 


3.72 


-0.07 


-0.16 


4.19 


254 


0.27 


3.99 


-0.07 


-0.17 


4.45 


165 


Pd2MnIn 


0.11 


3.81 


-0.05 


-0.14 


3.98 


279 


0.17 


4.10 


-0.05 


-0.15 


4.26 


173 


Ni2MnSn 


0.21 


3.74 


-0.06 


-0.13 


4.10 


321 


0.21 


4.00 


-0.06 


-0.14 


4.35 


251 


Pd2MnSn 


0.09 


3.79 


-0.06 


-0.13 


3.91 


368 


0.07 


4.07 


-0.06 


-0.14 


4.15 


263 


Ni2MnSb 


0.15 


3.76 


-0.03 


-0.09 


4.04 


309 


0.22 


4.04 


-0.02 


-0.06 


4.46 


234 


Pd2MnSb 


0.08 


3.82 


-0.03 


-0.08 


3.95 


193 


0.09 


4.13 


-0.03 


-0.07 


4.29 


296 


Ni2MnTe 


0.22 


3.83 


0.03 


0.01 


4.31 


120 


0.21 


4.04 


0.02 


-0.03 


4.48 


46 


Pd2MnTe 


0.14 


3.89 


0.04 


0.04 


4.21 


-172 


0.13 


4.17 


0.03 


0.02 


4.45 


-154 



We choose two different lattice parameters for each sys- 
tem: 6.07^1 and 6.38A. The first corresponds to the lat- 
tice constant of Ni-based alloys while the second is the 
characteristic for the Pd-based systems. The both sys- 
tems are experimentally well studied. In Fig.[5]we present 
spin projected total DOS and sp-electron DOS of the X 
and Z atoms for the lattice constant of 6.38A. Note that 
Ni {3d) is isoelectronic to Pd (4d). The calculated mag- 



netic moments are given in Table HT] The comparison of 
the results obtained for the same lattice constant reveals 
similarity of some features: the value of the Mn mag- 
netic moment, sp-electron spin polarization and DOS of 
the sp states of X and Z atoms around the Fermi level 
(see Fig. [8] and Table HI)) . On the other hand, well below 
the Fermi energy the DOS of the Ni-based and Pd-based 
systems differ strongly. The peaks of the DOS of the Ni- 
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based compounds are higher and narrower than in the 
Pd-based systems. This is mainly the result of a more 
delocalized character of the Pd Ad orbitals compared to 
the Ni 3c? orbitals. As we will show, despite the strong 
difference of the occupied part of the DOS the magnetic 
properties of two systems are rather close that gives ad- 
ditional support to the conclusion that the states lying 
close to the Fermi level play the most important role in 
the formation of magnetic properties. 

In Fig. [8] (left panel) we present the parameters of the 
Mn-Mn exchange interactions for Ni2MnZ and Pd2MnZ 
(Z=In, Sn, Sb, Te) for two different values of the lat- 
tice constant. The Curie temperatures estimated within 
mean-field approximation are given in Table |TT1 For 
both systems the patterns of exchange parameters for 
the equal lattice constants are rather similar. The pa- 
rameters are very close to each other for Z = Sn and 
Sb. For Z = Te, the difference is stronger. This differ- 
ence is more pronounced for near neighbors. For example 
the difference in the first two exchange parameters makes 
Pd2MnTe antiferromagnetic in contrast to the ferromag- 
netic Ni2MnTe. 

The analysis in terms of the competition between two 
types of exchange interactions suggested in the previous 
section is helpful also in this case. The long range behav- 
ior of the exchange parameters, specifically the RKKY- 
type oscillations can be related to sp-electron spin po- 
larization and total sp-electron DOS at the Fermi level. 
As seen from Table |lT] and Fig. [5] these quantities are 
very close to each other in both systems that explains 
the similarity in the long range behavior of the exchange 
interactions. However, the short range behavior depends 
very much on the contribution of the antiferromagnetic 
superexchange. To reveal the importance of this mech- 
anism we present in Fig. [5] the total DOS (spin-up plus 
spin-down) above the Fermi energy for both systems and 
for the lattice constant of 6.38A. As seen in Fig. [51 the 
DOS above the Fermi level is not identical for X = Ni 
and X = Pd. Therefore the deviations in exchange pa- 
rameters at least within few coordination spheres can be 
expected. A very strong similarity in the patterns of the 
exchange parameters for Sn and Sb alloys compared to 
the In and Te systems cannot be explained on the basis 
of the differences in the DOS and might be accidental 
resulting from the cancellation of the contributions of 
different states to the superexchange. 

Next, we comment on the volume dependence of the 
magnetic properties. As seen in Table HIl and Fig. [51 the 
compression of the volume (the reduction of the lattice 
parameter) leads to the decrease of the Mn magnetic mo- 
ments and an increase of the absolute value of the lead- 
ing exchange parameters. The reduction of the magnetic 
moments is an expected result which is a consequence 
of the increased interatomic hybridization and broaden- 
ing of the electron bands. The volume dependence of 
the exchange interactions is less straightforward. In gen- 
eral this dependence is non-monotonous and in concrete 
physical situations one can find both an increase and a 
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FIG. 9: (Color online) The total DOS of Ni2MnZ and 
Pd2MnZ (Z=In, Sn, Sb, Te) above Fermi level. 



decrease of the exchange interactions with volume con- 
traction. A detailed discussion of these aspects can be 
found in Ref. ]48| where the electronic structure, exchange 
interactions and Curie temperature of the fuU-Heusler al- 
loy Ni2MnSn are studied as a function of pressure. It was 
shown that in low pressure region the Tc behavior is in 
qualitative correlation with empirical interaction curve 
of Kanomata et. al. which describes the dependence of 
the Curie temperature of the Mn-based Heusler alloys on 
the Mn-Mn distance»2^ In agreement with experiment we 
have found that at ambient pressure Tc increases with 
increasing pressure that is dTc/dP > 0| 95i96,97 rpj^^ pres- 
sure dependence has a maximum at 3.6 A. Indeed, as 
seen in Table [IT] in agreement with experiments we ob- 
tain the same behavior in the Tc (Tn for Pd2MnTe) for 
all compounds except Pd2MnSb in which the Curie tem- 
perature decreases with reduction of Mn-Mn distance. 
To summarize, the magnetism of the Heusler alloys con- 
taining different d {3d or Ad) electrons appeared to be 
qualitatively similar whereas quantitative differences re- 
sult from the stronger delocalization of the Ad electrons 
and larger lattice parameters of the compounds having 
Ad elements. 



VI. SUMMARY AND CONCLUSIONS 

We performed a systematic first-principles study to re- 
veal the exchange mechanisms in various Mn-based semi- 
and full- Heusler alloys. The calculation of the exchange 
parameters is based on the frozen-magnon approach and 
the Curie temperature is estimated within the mean-field 
approximation. Due to large separation of the Mn atoms 
and the local moment nature of magnetism in these sys- 
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terns the exchange couphng is indirect and is mediated 
by the conduction electrons. The results obtained are 
interpreted using s-d mixing model of Anderson. To un- 
derstand the dependence of physical characteristics on 
the valence electron number we go beyond the stoichio- 
metric compositions employing virtual crystal approxi- 
mation. The influence of non-magnetic 3d versus 4c? elec- 
trons on exchange coupling is discussed. 

We show that magnetism in these systems strongly de- 
pends on the number of conduction electrons, their spin 
polarization and the position of the unoccupied Mn 3d 
states with respect to Fermi level. Various magnetic 
phases are obtained depending on the combination of 
these characteristics. The magnetic phase diagram is 
determined at T = 0. We find that in the case of a 
large conduction electron spin polarization and the un- 
occupied Mn 3d states lying far above the Fermi level, 
an RKKY-type ferromagnetic interaction is dominating. 



On the other hand, the antiferromagnetic superexchange 
becomes important in the presence of large peaks of the 
unoccupied Mn 3d states lying close to the Fermi en- 
ergy. The resulting magnetic behavior depends on the 
competition of these two exchange mechanisms. The ob- 
tained results are in good correlation with the conclu- 
sions made on the basis of the Anderson s-d model and 
with available experimental data. These findings suggest 
that a targeted influence on the corresponding physical 
quantities can provide a useful tool for the fabrication of 
materials with desired physical properties. 
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